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We study curvature perturbations in the anisotropic inflationary model with a complex scalar 
field charged under a f/(l) gauge field in Bianchi I universe. Due to Abelian Higgs mechanism, the 
gauge field receives an additional longitudinal mode. We verify that the dominant contributions into 
statistical anisotropies come from matter fields perturbations and one can neglect the contributions 
from the metric perturbations. It is shown that the contribution of longitudinal mode into the 
statistical anisotropy power spectrum, though exponentially small, has an opposite sign compared 
to the corresponding contribution from the transverse mode. 

PACS numbers: 



I. INTRODUCTION 



Simple models of inflation predict almost scale invariant, almost adiabatic and almost Gaussian perturbations on 
Cosmic Microwave Background (CMB) which are in very good agreements with cosmological observations 1 . There 
may be indications of statistical anisotropies on CMB [5J |3] which can not be generated in simple models inflation 
based on scalar flelds. Although the statistical significance of the possible statistical anisotropies on CMB is not high, 
nonetheless this opens up the interesting possibilities that primordial seeds in generating curvature perturbations 
during inflation may not be statistically isotropic. This will shed new light on the mechanisms of inflation. 

One can parameterize the statistical anisotropy via [1] P({k) = Po{k) (l + f/* cos^ 6) in which P(;{k) represents the 
curvature perturbations and 9 is the angle between the preferred direction in the sky which breaks the rotational 
invariance and the momentum vector k. Constraints from CMB and large scale structure indicate that \g^, \ < 0.4 

Motivated by these observations, there have been many attempts in the literature to generate primordial anisotropies 
during inflation. The natural way to break the statistical anisotropy during inflation is to employ a gauge field or a 
vector field to seed the anisotropies at the order of few percent which may be detectable on CMB [7H27] . 

As demonstrated in |28j . in models employing a massive vector field in which the gauge invariance is broken 
explicitly, the excitations contain a ghost which is not acceptable physically. Therefore, it is crucial that the vector 
field is protected by a gauge symmetry so the longitudinal mode of the vector field excitations is not physical. On 
the other hand, because of the conformal invariance of models with gauge fields, any excitation of gauge field during 
inflation is diluted and can not seed the desired anisotropies. Therefor it is essential that one breaks the conformal 
invariance while keeping the gauge symmetry explicit. This approach was employed in different contexts in |29H46| . 
Specifically, one can consider the model in which the U{1) gauge kinetic coupling is a function of the inflaton field 

with the action AC = "■^^'^^ F^i^F^^'^ in which (f) is the inflaton field and i^^^ is the U(l) gauge field strength. If 
one choses /((/>) such that /((/)) cx a^^ then a constant source of electric field energy density is turned on at the 
background level and the gauge field quantum fluctuations remain scale invariant. As shown in [3S] the inflationary 
system admits an attractor solution in which the anisotropy reaches a small but cosmologically detectable level. 
Cosmological perturbation for this model in which inflaton field is a real scalar field with no charge coupling to the 
gauge field is studied in great details in [47H52] . 

In this work we perform the cosmological perturbation theory for the model presented in |39| in which the inflaton 
field is a complex scalar field charged under the U{1) gauge field with the electric charge coupling e. Under the 
Abelian Higgs mechanism, the gauge symmetry is spontaneously broken and the gauge field acquires a dynamical 
mass in the form e^p^A^A^^ in which p is the radial component of the complex inflaton fleld </>. As we shall see this has 
interesting implications for the cosmological perturbations and in generating statistical anisotropies. Namely, as in 
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usual Abelian Higgs mechanism, one scalar degrees of freedom is eaten by the gauge field and the longitudinal mode 
of the gauge field is excited. As a result, along with the two transverse modes of the gauge field, the longitudinal 
excitations will also contribute into anisotropy analysis. 

The rest of paper is organized as follows. In Section ITT] we present our model and classify the metric and matter 



perturbations. In Section III we present the second order action which will be used to calculate the anisotropic 
curvature perturbations power spectrum in Section |IV[ In Section |V] we demonstrate that the leading perturbations 
come from the matter sector. Summary and conclusions are given in Section |VI[ We relegate technical discussions 
about the choice of our gauge, integrating out non-dynamical fields and the detailed forms of the second order action 
into Appendices. 



II. ANISOTROPIC INFLATION FROM CHARGED SCALAR FIELD 



Here we present our model and the metric perturbations and gauge choice in Bianchi I background. 
The model we are interested in is studied at the background level in |39j . It contains a complex inflaton field 
which is charged under the U{1) gauge field with the electric charge e. The action is 



S 



A/2 



2 2 " 4 

in which Alp is the reduced Planck mass and the covariant derivative is given by 

D, 



As usual, the gauge field strength is given by 



(1) 



(2) 



(3) 



As explained above, we have inserted the time-dependent gauge kinetic coupling in order to break the conformal 
invariance such that the gauge field excitations acquire a nearly scale-invariant power spectrum and survive the 
exponential expansion. In order to obtain a scale- invariant gauge field power spectrum one requires / oc a"^. This 
corresponds to a constant electric field energy density during inflation. 

We assume the model is axially symmetric in field space so V and /(</>) are only functions of 4>4) — It is more 
instructive to decompose the inflaton field into the radial and angular parts 



cl>{x)^p{x)e^'^^\ 

soy = l/(p)and/2(<^) = /2(p). 

As usual, the action ([l]) is invariant under local gauge transformation 



Af_i — > A^ 



^df,e{x) , 



e{x). 



In terms of p and 9 the action ([T]) is given by 

5 = 



Pip) 



R - ^d^.pd'^P ~ y [d^e + bA^) {d^e + e^'') - ^F^.F^-- - V{p) 



(4) 



(5) 



(6) 



We assume that the gauge field has a non-zero classical value along the x-direction so A^ — (0, 0, 0). As a 

result, the background space-time is in the form of type I Bianchi Universe with the metric 



ds^ = -df + e2"(*) (e-4-(*)dx2 + e^^^'\dy^ + dz^)) 



(7) 



In this view a measures the averaged Hubble expansion while a{t) measures the level of anisotropy. In order to be 
consistent with cosmological observations, the level of anisotropics should be very small so a/a 1. 
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h4(T 



A 



The background fields equations are given in [35] 

a*(/'(p)e" 

p + iap + Vp+ {-f{p)f,p{p)A^ + e'pA^) e-2"+4- 
1 , fl 



= 



V{p) 
V{p)- 



-^f{p)A' + ^A' 

lf{p)A' + ^A' 
3 

lf{p)A' - ^A' 



(8) 
(9) 
(10) 

(11) 
(12) 



in which a dot indicates derivative with respect to t. 

As in conventional models of inflation, the background expansion is driven mainly by the potential term V. However, 
the gauge field also contributes in the background expansion in the form of electric field energy density turned on 
along the x-direction. In order for the anisotropy to be small we require that the electric field energy density to be 
very small compared to V. This corresponds to i? <C 1 in which 



A^fipfe^ 
2V 



(13) 



The attractor solution 



It is more convenient to express the background metric Q in the following form 

ds^ = a(7/)2(dr;2 + dx^) + b{r^f{dy^ + dz^) 



(14) 



in which a 
slow-roll parameters 



and b 



Here we have defined the conformal time dr] via dt = a{rj)dri. Let us define the 



H 



H 

2HH 



en = 2HeH{2tH - Vh) 



(15) 



We are working in the slow-roll limit in which eH-^TlH ^1- To leading order in slow-roll parameter and anisotropy 
a~6~-l/iJ77. 



Although the anisotropy is very small, ^ 1, so the Hubble expansion rate in modified Friedmann equation ( 10 ) is 



mainly dominated by the isotropic potential term, but the back-reactions of the gauge field on the inflaton field induce 
an effective mass for the infiaton as given by the last two terms in Eq. ([9]). This in turn will affect the dynamics of 
the inflaton field. As shown in [351 the system reaches an attractor solution in which R (x en- For this to happen we 
need f{p) oc a" with n ~ —2. Indeed, the background expansion is given by 



So if one chooses 



a (X exp 



/ c>c exp 



dp 



V 



MlVp 



dp 



V 
MlVp 



(16) 



(17) 



this yields f oc a"^. 

The exact form of / therefore depends on V{p). For the chaotic potential used in 



we have 



V 



1 



1 2 

m p 



fip) = exp 



2Ml 



(18) 



with c a constant very close to unity. As shown in '38' during the attractor phase the effective inflaton mass is reduced 
by the factor factor 1/c such that during the attractor phase the inflaton evolution is given by dp /da ~ —MpV^p/cV. 
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The cosmological perturbations for this background was studied in details in [47H52] with the conclusion that in order 
not to produce too much anisotropy one needs c — 1 ~ 10~^. 

For our model, following j39| . we consider the symmetry breaking potential which is physically well- motivated for 
the charged scalar field in the light of Abelian Higgs mechanism. The potential is 

in which A is a dimensionless coupling. The potential has global minima at fi — ±Af/\/A- The inflaton field rolls near 
the top of the potential so in the slow-roll limit, the potential can be approximated by 

V ~ ^M^ - — p2 _ (20) 
4 2'^ V ; 



From Eq. ( 16 ) we have 

M 

'p 



aocp-/^ , P.^J^- (21) 



To have a long enough period of slow-roll inflation we require pc 3> 1. 



Motivated by this, from Eq. (171 we see that to find an attractor solution with a near scale invariant gauge field 
power spectrum (i.e. a scale invariant electric field power spectrum) we take |39j /(p) oc p^^ with p very close to Pc- 
Noting that p cx a^^^P" cx (—77)^^/^= this yields 



2c 

V \ _ P 



in which r]e indicates the time of end of inflation. We assume that at the end of inflation / reaches its canonical value 
/(te) = 1 and the isotropic FRW universes emerges at the end of inflation. As we shall see, the strength of anisotropy 
is measured by the parameter / given by 

I^'—1 = P^. (23) 



During the attractor phase 



i?^^ , (24) 
2 ' 3 3 ^ ^ 



This indicates that the anisotropy is at the order of slow-roll parameter during the attractor phase. 

In this picture inflation ends when the back-reaction of the gauge field on the inflaton field via the interaction 
e^p^ApA'^ induces a large mass for the inflaton. Comparing this with the inflaton mass M, inflation ends when 
e^e~^"= A^(?7e) ^ in which ae indicates the number of e-folds at the end of inflation. As shown in [39] the end of 
inflation depends logarithmically on e. More specifically, noting that during the attractor phase [HI cx e^", we 
obtain 

Ine , ^ 

ae + ... (25) 

where dots indicate the dependence on other parameters such as Pc and the initial value of the gauge field. As one 
expects, the larger is the gauge coupling e, the shorter is the period of inflation. This is easily understood from the 
induced mass term e^A^A'^p^ for the inflaton field due to Higgs mechanism. 



B. Perturbations 



Now we look at the perturbations of the background metric ( 14 ) . Because the gauge field has a component along the 
x-direction, the three-dimensional rotation invariance is broken into a subset of two-dimensional rotation invariance 
in y — z plane. Therefore, to classify our perturbations, we can look at the transformation properties of the physical 
fields under the rotation in y — z plane. As mentioned in [17W33] the metric and matter perturbations are divided 
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into scalar and vector perturbations for a general rotation in y — z plane. It is also important to note that there are 
no tensor excitations in two dimensions. 

The most general form of metric perturbations is 



( 



\ 



'2c? A a^d^jS ab{d,B + Bi) \ 



(26) 



Here A, f3, B,ip,j,Tp and E are scalar perturbations and Bi,Ti and Ei are vector perturbations subject to transverse 
conditions 



0. 



(27) 



In Appendix |A] we have presented the properties of metric perturbations under a general coordinate transformation. 
In our analysis below we chose the following gauge 



(28) 



which from Appendix [A| one can check that it is a consistent gauge. Note that the gauge (281 is similar to the flat 
gauge in standard FRW background. 

As for the matter sector we choose the unitary gauge 9 = Oso (j)is real. Also, exploiting the two-dimensional rotation 
symmetry we choose A;z = so k = [k^, ky, 0) in Fourier space. Therefore the scalar and vector perturbations of the 
matter sector, 6A^^^ and are 



Mjf) = {5Ao,5A^,dyM,0) 



54^) = (0,0,0, i^). 



(29) 



With these decompositions of the metric and matter fields into the scalar and vector sectors, one can check that 
these modes do not mix with each other and one can look at their excitations and propagation separately. In this 
work we concentrate on the anisotropics generated from scalar excitations which are more dominant compared to the 
anisotropies generated by vector excitations. Therefore, for the rest of analysis we set D = Ti = Bi = Q. 

III. SECOND ORDER ACTION 

Here we present the second order action for the scalar perturbations. Our is to find the second order action both 
for the free fields and for the interactions. As we shall see the fields 5Aq,A^P and B are non-dynamical in the sense 
that they have no time-derivatives in the action. As a result, their equations of motion give constraints which can be 
used to eliminate them in terms of the remaining dynamical fields Jp, 5Ai^ M and 7. 

The second order action for the scalar perturbations is 



drjd X 



2hb'A,^f3^^ + ab{"- + ^)A^yB, y + abj^^yA^^y - a%'^V{pn)A^ - yfo^o^^^^ - '^P,xyB,^y 



ab 6^ 6^ b' a' b"^ 6^ 



i2 a2 2 2 
-,2 



2a2 



{Kflly 



b\b' 



(— - —)Xxy + -w^P - b P'oASp' - b'pof3,Jp,, - abp'oBJp^y 



4 ' b 



-Sp'. 



a 



-5pl 



b'ptSAf, - e'b'piA^p^JAa 



+e^abplAx'y,xyM^y 
pb^ 



o^pImI - e^b'^p^A^ASAi - e^b^poAlASp 



-M'iMo,x 



~ '2 

f2f^2 ^-252 

^A'^ASA\ + ^A'^ASAo 



pb^ 
2a2 



^-y'ptSAi - ^b'AiSp' - 2e%^poAx6p6A, 



SA'l 



2a2 -^^o.- 



Kl,xyM[y + '—A'^^,,y5AQ 



A'B,ySAi 



'—A'^B^yM^Xy 



f2 

2 y 



^ - pM'JA^^y - L.dAly ~ i^Mly 



fSAl^yM^Xy 



a2 



Jf^AiJA[Sp 



-2'-^A'JAo,Jp^ 



—^ASP 



^'"^^'^ 6 p^ 



2a2 



f f,Ppb^ j^2^p2 _^^b'^ X„2 „2;,2i 



2a2 



-V^ppSp^ - a^b'^VpASp 



(30) 
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in which a prime indicates derivative with respect to conformal time. 

As mentioned above, the excitations SAQ,A,f5 and B have no time-derivatives so they are non-dynamicaL The 
details of ehminating the non-dynamical excitations in terms of dynamical perturbations are given in Appendix [B| 

The final second order action is a complicated function of Sp, SAi , AI a nd S pecifically, integrating out SAQ,A,f3 



and B one encounters the functions and as defined in Eqs. (B17) - (B38|. At this level it seems hopeless to 
get any insight into the form of the action and the prospects for analytical analysis. Happily, the analysis becomes 
considerably simple if one notice the following effects. Looking at the formulae for A^ and A^ it is evident that Ai is 
the key parameter which controls the form of other A; and A^. Now let us look at the function Ai 

Following the procedures of integrating out the non-dynamical fields in Appendix |b] one can check that Ai comes from 
integrating out SAq. Neglecting the anisotropy for the moment, the ratio of the second term in Ai compared to the 
first term scales like e^a^//^ ~ r]^ / rf . Therefore, during the early stages of inflation in which —77 3> —r]e, the 
second term in Ai is completely negligible compared to the first term. In this limits all Aj and A^ collapse to simple 
forms and we will be in the limit somewhat similar to [49] . In this limit the effect of gauge coupling e is sub-dominant 
in the action and, as in ^49) , the mode excitations are massless. On the other hand, as inflation proceeds the second 
term in Ai eventually dominates and we enter the second phase in which all Ai and Ai are proportional to e^. The 
time when the two terms in Ai become comparable, denoted by -q^^ is given by 

1/3 



In order to obtain analytical results, we solve the mode function in each period, for r] < rjc and for rjc < rj < rje and 
glue the solutions by appropriate matching conditions. 

Having this said, one may wonder why Ai plays such a prominent role. The answer to this question is provided in 
Section|V]in which we demonstrate that the leading interactions come from matter sector. So it is not surprising that 
only Ai, which originates from integrating out 6Aq, will have a prominent efect while the other parameters A^ and Ai, 
which have their origins in integrating out metric flclds A, j3 and B, are negligible. 

A. Second Order Action in the Slow roll approximation 

Here we obtain the second order action in the flrst order of slow-roll and anisotropy parameter /. Including the 
first slow-roll and anisotropy corrections into the background expansion one can check that 

a ~ H-\-T^)-^-'" , bc^H-^{-ij)-^~'-^' (33) 

Also, as mentioned above, we chooses the coordinate system such that 

k = {kx,ky,0) , kj: = k cos 9 , ky = -ksinO . (34) 



From the above formulae, and using Eq. (22) for the function /(ry), one can obtain the following expressions which 
would be useful later on 

" (-ry)-i(l + e^) ^^(_,^)-i(l + e^ + /£^) 



a 



b 



° (-77)-2(2 + 3eff) y = (-77)-2(2 + 3eH + 3/eH) 



a 



J - {-vr\- sin' OIch) ^ - i-vr'{~ sin' Olen) 

^ = {-Vy\-2-2eH-VH + 2IeH) ^ = (-7,)-'(2 + 96^ - 37?^ + 6/e^/) . (35) 

So far our analysis were generic and we have not used any particular form of the potential. From now on, we 



concentrate on the symmetry breaking potential given in Eqs. (19) and (20) and consider the phases rj < rjc and 
rjc < T] < rje separately. 
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First we consider the period in which rj < rjc so the term containing e in and Xi are neghgible. Considering the 
leading corrections from the slow-roll and anisotropy expansion yields (for details see Appendix O 



82^'' — J dr] (Pk^Lpp + L^j + Ld^Di + LD2D2 + I^pi + -^pDi + LpD2 + Lj^-^ + L^^^ + i/jj 



in which 



L 



77 



12- 



:(1 - 2sin2 



I- I I- I' 
+ {~vy^ ( 2 + 15€h + leni-Q + H sin^ 9 + cot^ ( 



1\ 62 



1 - / 



1 



e/sin^ 6l/(^(5p*i:ii +C.C.) 



cos^ 61 £'2 + C.C.) 



^ k^f^ilen sin'*6lcos2 6' ) (7*1)1-^0.0.) + 
63 1 



1 

2a2 71 



7 I?i + c.c 

1\ / M2 



2 k^ f ^/SR^ sin^ 9 cos^ e 



1 -I- c.c. 



, , -k^f^ySIeH sin^ 61 cos^ 6'(1 -f cos^ 61) ) (-f'*D 
4a^ 77 /V 

fcVySie^sin'' ecos^ 61^ (^7'*i:)i 4- c.c.) , 

2J-2cos2(? 



ry / \ 16a/ 
\/263 



cos f sm 



/2) 



(3 -I- cos2 



7 *L'2 + c.c. 



4a 



A:2ve]Fsin2 6'cos2 6'j ((5p'*(57' c.c.) -f A;2(5p* 57 c.c.)) 



Here Di and _Di stands for the transverse and longitudinal modes respectively defined via 

Di = SAi - ik cos 9M 

D2 = cos 96 Ai + ik sin^ 9M . 

Furthermore, we have defined the canonically normalized fields 



a kie 



Ik 
D^k 
'Wk 



X y 

Ik = 



a k^ky 
V^~k^ 



Zk 



V2 fc2 
-fsin9Dik = -fs\n9vk 



2V2Xk y l-I 



EH _ eM 

:bD2k = 



2V2Xk y l-I 



bwk 



(36) 

(37) 
(38) 
(39) 
(40) 
(41) 
(42) 



(43) 
(44) 

(45) 

(46) 
(47) 

(48) 
(49) 

(50) 

(51) 



Note that Eqs. (37)-(40) represent the free-field actions for Sp,j,Di and D2. As expected, during this phase in which 
the effect of gauge coupling e is sub-leading, similar to [49 , Eqs. (37)- (40 1 represent nearly massless fields with almost 
scale-invariant power spectrum. 

The interaction terms are given by Eqs. (41 )-(45 ). For technical reasons the interaction terms are presented in terms 
of the original non-canonical fields. At first it seems that Eq. (45) is the leading interaction compared to the other 
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interaction terms. However since 7 = ^^^5^7, this term does not have any 6'-dependence so it only re-normahzes 
the isotropic curvature perturbation power spectrum. Furthermore, by using the In- In formahsm in Section |IV| one 



can easily check that its contribution is suppressed compare to Eq. (41 ) . The reason is that this term is proportional 
to {6p *6j +C.C.) and when we write it in terms of the canonical fields the leading terms cancel each other. Therefore 
we do not study the contribution from this term here. Also note that during this phase of inflation the interaction 
LpD2 between the longitudinal mode and the infiaton field is much smaller than the interaction Lp^-^ between the 
infiaton field and the transverse mode. As a result the most important interaction term during this period is given 
by Eq.(41) from the transverse mode Di. 



The next step is to calculate the mode functions related to our canonical variables (5p, 7, -D2- As mentioned 
above, since during this phase the modes are nearly massless, we have the following standard form of mode functions 
with the Bunch-Davics initial condition 

m.-k^^f^^'^^l-f )• (52) 



As we mentioned below Eq. (31) during the period rjc < t] < rje the effect of the ga uge c oupling e will become 
important. During this phase the dominant contributions in and Aj in Eqs. (B17) - (B38l come from the terms 



containing e. Expanding to leading order in terms of the slow-roll parameters and / the second order action is 



d-qd^k Lpp + LdiDi + L 



Where, 



^pp 



L 



DiDi 



L 



D2D2 



-ID,' 

1 ' 



D2D2 



4A 



4A 



1 

7V 
1 

1 

p^2 



L 



-'pD2 



a 
a 



Ve/sin^ ef[5p*D[+c.c 
VgI cose f(6p*D'2 + c.c 



During the second phase the canonical variables are defined as 



ic (m) 

-sm9fD,k^- sine fvi'''^ 
a a 



/■J-, ^ r (m) 

fD2k = -fwl 



L 



pD2 



\Sp\' 



\D,\' 



(53) 

(54) 
(55) 
(56) 
(57) 
(58) 

(59) 
(60) 

(61) 



The crucial difference compared to the first phase is that the modes are not massless anymore. Due to interaction 
e^p^Ap_A^ the scalar perturbation Sp, the transverse perturbati on Di an d the longitudinal perturbation D2 acquire 
time-dependent masses proportional to as can be seen in Eqs. ( 54 )-( 56 ). Unlike the previous phase, the longitudinal 



mode also plays important roles as we shall see below. Note that in above equations the superscript m in u^™^ , w^™'' 

and if^™'' indicates that these modes are massive. 

As in the first phase, for the purpose of the calculating the curvature perturbations power spectrum, we look into 
interactions between Sp and other fields. As before, the interaction Lpj does not have any directional dependence so 
we have not considered it in above action. Therefore we are left with Lpjj-^ and Lp^^ ■ 
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Comparing Eg. (48) with Eq. (59), we see that, although Sp — Sp and Di — Di, however D2 is not equal to 
Z?2- Also Eq. (57) is proportional to sin^ 6* while Eq. (58) is proportional to cos 6*. As a result we can guess that 



the contributions of the longitudinal mode in has a different sign than the corresponding contributions from the 
transverse mode. So the question arises whether or not we can produce a positive factor from the longitudinal 
mode (from [49] we know that g^, is negative for the transverse modes). We will come back to this question when we 
calculate the power spectrum of curvature perturbations. 

As in the first phase, the next step is to calculate the mode functions. The Sp^^ mode function is. 



"(-k)«p(-k) 



CikH, 




C2kH. 




n : 



M4 



-Ve- 



(62) 



Here H^^^fix) are the Hankel functions with rank | and Cik are constant of integrations to be found by matching 
conditions at rj — rj^. 



Also Z?ik mode function is 



■ sm ( 



-T] 



bsmOf 



(m). 



(m)* , f 
(-k) p(-k) 



dikH. 





A 



(63) 



Finally the Z32k mode function is 



f 1 / (™) I (™)* t 

/) [^k cpk + w;(_k)c;(_k) 



eikH. 



(1) 
3/4 



277 



2 ) + ^2^^/3/4 [ ^ 



(64) 



Our goal is to find the constants of integration Cik,dik,eik , i = 1,2 with the appropriate matching conditions at 
rj — rjc- To impose the matching conditions we require that the original fields 6p, Af, Ai and their derivatives to be 
continuous at rj = rjc- Imposing the matching conditions one obtains 



Cl,2 



dl.2 = 



ei,2 



TT 


e 


-ikric 


8i 


V- 


-2kr]c 


TT 


e" 


-ikrjc 


8i 


V- 


-2kr}c 


TT 


e" 





±H. 



(2,1) 
3/4 



±H. 



(2,1) 
3/4 



8i ^J—2kric 



(2,1) 
3/4 





(65) 
(66) 
(67) 



IV. POWER SPECTRUM OF THE CURVATURE PERTURBATIONS 



Having calculated the final outgoing mode functions, we are ready to calculate the curvature perturbation power 
spectrum. We are interested in anisotropics generated in curvature perturbation power spectrum. The anisotropics 
are generated by interactions Lpjj-^ and ipzjj from the coupling of the transverse and longitudinal mode to 5p. The 
corresponding Feynman diagrams are given in Fig. 1. 

Using the standard In-In formalism for the curvature perturbation power spectrum |54H56| we have 



< 5p ive) 



Texp 



Hi{f]')drj' 



dp (ry) 



Texp 



Hi{T^')dri' 



(68) 



where T and T respectively denote the time-ordered and anti-time-ordered products and Hi refers to the interaction 
part of the Hamiltonian in the interaction picture. As for 770 we can take 770 — > —00 so the modes of interests were 
originally deep inside the horizon. 
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5p 



D1 



5p 



5p D2 5p 



(a) 



(b) 



FIG. 1: The transfer vertices for the interactions of the inflaton field Sp witli the gauge field excitations Di and D2. The left 
figure represents Lpo^ as given by Eq. (571 while the right figure represents I/pOj given by Eq. (58 1. 



To leading order the contribution of anisotropy in inflaton power spectrum, A < 6p [rje) >, is 



A < (5/9 > = 



»)e r-rii 

drji / dri2 

rio J no 



(69) 



In our model the form of the interaction Hamiltonian in two phases are different. This is because during the first 
phase only the interaction between the transverse mode and inflaton, Lpu-^ , is important while during the second 
phase both interactions from the transverse and longitudinal modes, ipD^ and Lpjj^, are important. 
We can connect these two forms of Hamiltonian into each other by using the step function as follows 



l) (^^SJsin^e^ {6p*D[+c.c.) 
JIJ V a 



^'^VeismH^ {Sp*D[+c.c.) ~ 



6/ cos 6* {Sp*D'2 + c.c.) 



(70) 
(71) 

(72) 



Here Hf and refer to the Hamiltonian in the first and second phase respectively. The terms containing Di are 
the same in both Eqs. (71 1 and (72 1 while the terms containing D2, coming from the longitudinal mode, appears only 



in Eq. ( 72 ) during the second phase. 



By plugging back Eq. ( 70 1 into the Eq. ( 69 ) we have 
^ 9 



A < (5/9 > 

< iVe) > 



192/ 



rn 

drji / dr]2 

J Vo 



«(•)! )«■"'(>?) 



Im 



X 9{ric — '7i)^(?7c — V2) + sin^ 9Im 
X e{r]i - r]c)9{T]c - m) + sin'* 6Im 
X 6{rii — ric)0{r]2 — rj^) + cos^ 6Im 



M™(ryi)«"*(77) 
M™(ryi)«™*(77) 
7."(77i)u™*(r7) 



Im 
Im 
Im 



u{r]2)u"^* {ri)v* (771)1; (772) 



X 6l(r/i - Vc)0{r]2 - ?7c) , 



(73) 



where |u^''-'(7ye)|^ 



21^3^2 is the amplitude of the free inflaton field fluctuations. 
So the correction to the power spectrum has four terms as follows 



Pi = 



192/ 



P9 



24/sin2 0N^ 
192/ 



l-Vl 

drji / dri2 
J no 



U(77l)u™*(77) 



Im 



Im, 



u{r]2)u"^* {r])v* {r]i)v (772) 



7.(772)7^"* (77)7;"* (77i)« im) 



(74) 



24/sin2 e { 2NcNe ~ 2N^ ~ 



(75) 
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I — — 1 sm Olr 



2 N 
241 sin' e{N^ + Nl + -N, - 2N,N, - 



7/'"(77i)«™*(r7) 



Im 



(76) 



P4 



192/ 



dm / dTj2 I ^ ] cos^ eim 



Im 



27r(3/4)4 \^'^ 2 

As a result the fractional change in the curvature power spectrum due to anisotropy, to leading order, is 



(77) 



P1+P2+P3+ Pi 

N, 



24/ N, 



(78) 



6 J ■ 2Tr(3/4)'^ 

In this formula, stands for the total number of e- folds which we take to be 60 and is the number of e-fold from 



sm^ + L X U^r^'V^ ] (TVe - cos' . 



the time rj = rjc till end of inflation. Using rj ~ —1/aH and Eq. (32) we have 

27Ve 1 / /3e^\ 2N, 



(79) 



The last approximation is valid for typical parameter values such that the logarithmic correction in Eq. ( 79 1 is order 
unity. Taking to be 60 to solve the flatness and the horizon problem, we have Nc ^ 40. 

The first term in Eq. (781 come from the transverse mode which is similar to [49'. However, the last term in Eq. 



( 78 1 is due to longitudinal mode which does not exist in models in which the inflaton field is real. This is a direct 



consequence of Higgs mechanism in which the gauge field has an additional longitudinal excitation. We note that the 
contribution of the longitudinal mode in the 5* factor is positive which is an interesting observation. However, since 
krj^ ~ Hrjf. ^ e^^°, the last term in Eq. (78) is exponentially small so the dominant contribution in comes from 



the transverse mode. 

Since sin^ 9 = 1 — cos' 9, the leading correction to anisotropy power spectrum in Eq. 

'241 N' 



(78) is 



(80) 



g* 

Now demanding that \g^ \ < 0.3, we find that / ~ 10^'"'. 

Interestingly our leading expression for in Eq. (80) is the same as in 05]. This can be interpreted as follows. 
Because of the interaction e'p'Afj^A'^ the inflaton fleld, the transverse mode and the longitudinal mode become 



massive. This results in the complicated Hankel function for their mode functions u 



(m) 



and 



However, one 



can check that the argument of these Hankel functions is smaller than unity until towards the end of inflation. As a 
result, all these modes stay nearly massless as in [33]. However, unlike [49J, we also have the additional contribution 
from the longitudinal mode. The interaction term responsible for the contribution of the longitudinal mode is given 
by the second term in which is turned on only during r/ > ry^- However, the initial amplitude of D2 is much 



smaller than Di as can be seen from the coefficients in matching conditions Eq. ( 66 1 and ( 67 ) . This is because the 
wave function normalization of D2 in the first phase originates from the sub-leading interaction e'p'Afj^A^ while that 
of Di comes from the leading term f'Ff^^F^'^. As a result, the contribution of the longitudinal mode in anisotropy 
power spectrum is hugely suppressed as compared to that of the transverse mode and up to leading order we recover 
the result of gS]- 

As observed in [5 2) the infra-red (IR) modes of the vector field perturbations remain frozen on super-horizon scales 
which accumulate to renormalize the background gauge field. As a result, this can lead to a large value of 5* unless 
one takes N„ ~ 60 as we have assumed here. 



V. THE ORIGIN OF THE LEADING INTERACTIONS TERMS 



Having calculated the anisotropic power spectrum through complicated procedure of integrating out the non- 
dynamical fields and approximating Ai and other A; and Ai, one may wonder what the origins of the leading interaction 
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terms ipUi and Lpjj^, or alternatively LpA^ LpM and La^Mi are. Are they coming from the metric perturbations or 
from the matter sector? Here we show that the leading interactions actually come from the matter perturbations. 

To show this first we integrate out SAq and then read off the interaction terms containing the matter perturbations. 
The leading terms in the matter sector coming from integrating out 5A^ are 



A, 2a3 , + 

hi 1.3 

— sin^ cos 9(iM'6A* + c.c.) 
4a'* Ai 



k'^ cos sin^ 



Ai 



2a3 7/ 



f{i6p*M' + c.c.) 



(81) 



On the other hand, the leading terms for the matter perturbations present in the original action (without integrating 
out any fields) are 



6/ 



a 7] 



f(Sp*SA[ + c.c.) + —^k^f e cosOUMSAl + c.c.) . 



(82) 



So by adding Eq. (82) and (81 1 wc can obtain all the leading interaction terms for LpA-^ , ipM and La-^m as, 

f (tSp*M' + c.c. 



/6I J. k^cos^e b^ v67^3 
a rj Ai 2a^ rj 



k^fsin^9cos9] [ iMSA^ 



2a2 




Aa'i Ai 



2a^ 1] 

sin^ 61 cos 61 ) ( iM'6Ai + c.c 



(83) 



Interestingly, this is the whole leading action which was used in previous sections to calculate the anisotropic power 
spectrum. To see this more specifically, consider the phases r] < rj^ and r] > rjc separately. For the first phase, rj < rjc, 
Eq. (pi) yields 



^lead. 




6p*SA[ + C.c. 



k cos sin^ - V6If ( iSp* M' + c.c. 



k-^f^ sin^ 0COS0] iMSAl + c.c. 



2a2 



kf sin^ 61 cos 61 iM'6A[* + c.c 



Sp*D, 



c.c. 



(84) 



Interestingly, this is exactly the leading term interaction as obtained in Eq. (41 1. Similarly, for the second phase, 
7] > Tjc, Eq. (|83]) yields 



-^load. 




+ 




/ 6p*6A[ + c: 



fcV sin-" 6* cos 61 iMSAl + c.c. 



6/sin^ ef{Sp*D^ +C.C. + 



1\ 62 



6/ cos 



rj J a 



Of{Sp*D'^ 



+ c.c. 



(85) 



As expected, this expression is the sum of the leading interaction terms Eqs. (57 1 and (58). 

In summary we conclude that the leading interactions for generating anisotropics originate from the matter sector 
and one can neglect the metric perturbations in calculating the leading order corrections to the curvature perturbations 
power spectrum . Computationally, this is a very important result which considerably simplifies the perturbation 
analysis in similar models. This conclusion was also reached in !52]. 

This also explains why in t he processes of integrating out the non-dynamical fields only Ai plays prominent roles. 
As mentioned below Eq. (31) Ai originates from integrating out SAq which is the non-dynamical field in the matter 
sector. On the other hand, other A^ and A^ originate from integrating out the non-dynamical fields A, B and j3 in the 
metric side which should not play prominent roles as expected from the above results. 
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VI. SUMMARY AND DISCUSSIONS 



In this work we have studied cosmological perturbations in an anisotropic inflationary scenario with a complex 
scalar field charged under the U{1) gauge field. Because of the Abelian Higgs mechanism, the gauge field obtains 
the dynamical mass e^p^A^A''. As a result, the angular excitations of the complex scalar field is eaten by the gauge 
field so the longitudinal component of becomes excited. Furthermore, the interaction term e^p^A^Af^ induces a 
time-dependent mass for both the inflaton perturbations 5p and the gauge field excitations Di and D2 which become 
comparable to H at the end of inflation. 

The leading contributions to anisotropic power spectrum are given in Eq. (78). The first term comes from the 
interaction of 6p with the transverse mode, LpOi ■ In terms of Feynman diagrams this interaction is represented by 
the exchange vertex shown in Fig. 1 (a). This is similar to the result obtained in [49 . However, unlike j49], the 
longitudinal mode D2 also contributes into the anisotropic power spectrum. Interestingly its contribution into is 
positive as given by the last term in Eq. ( 78 ). In terms of the Feynman diagrams this interaction is represented by the 
exchange vertex shown in Fig. 1 (b). However, the longitudinal mode contributes only towards the end of inflation 
and its contribution to the anisotropic power spectrum is hugely suppressed compared to the contribution from the 
transverse mode. 

We also verifled that the leading interactions in anisotropic power spectrum come from the matter sector perturba- 
tions. In other words, to calculate the leading order corrections into the power spectrum, one can neglect the metric 
perturbations. Computationally, this knowledge simplifles the analysis considerably. This is particularly helpful when 
calculating the bispectrum and non-Gaussianities which we would like to come back in a future work. 

The issue of generating statistical anisotropy at the end of inflation via waterfall dynamics have been considered 
in |21j . However, it is shown in [22j that this mechanism does not work and the anisotropy produced purely from 
waterfall effect at the end of inflation is exponentially suppressed. This conclusion, however, was criticized in [23j . 
Having this said, we believe that the conclusion derived in [22] . which is obtained by a careful use of SN formalism, is 
valid. To get large enough statistical anisotropy in model of [H], one has to consider the evolution of the gauge field 
both at the background level and at the perturbation level during the entire infiationary period, as we did here. This 
point was also mentioned in |52| . We would like to pursue this issue in a future work considering the charged hybrid 
model using the standard In-In formalism as employed here. 
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Appendix A: Metric Perturbations 



Here we study the metric perturbations in Bianchi I background and their transformation properties under a general 
coordinate transformation. Consider the general coordinate transformation 



(Al) 



in which A and A are scalars and is vector subject to di£,\^=0. For the future reference note that by appropriate 
choice of A and A one can remove three scalar degrees of metric perturbations in Eq. (26) while the freedom from 
can remove only one vector degree of freedom. 
Under the coordinate transformation Eq. (Al) we have 



(0) 



(A2) 



in which ^^"^gafi is the background Bianchi metric given in Eq. (14). 
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More explicitly, one can check that 



A 


a ^ ' 


(A3) 






(A4) 


B 


a 


(A5) 




a 




1 


7 A - -A 

a 


(A7) 


V' 




(A8) 


E 


E-A 


(A9) 



and 



E, 



a 
b 



a 



(AlO) 

(All) 
(A12) 



Using the above transformation properties one can check that the following two scalar variables are gauge invariant 

-Tp (A13) 

(A14) 



1 



a ' 



Ha 



In this view (50^ represents the inflaton perturbations on -0 = surface which reduces to inflaton perturbations on 
flat slice in FRW background while 7 is identically zero in FRW background. 
In our analysis we adopt the following gauge 



iP^iP = E = E, = 



(A15) 



which one can check is a consistent gauge. Note that the three scalar conditions tp — ip = E = fixes three scalar 
freedoms A and A while th e vect or condition E^ = fixes the remaining one degree of freedom ^j^. The advantage 
in choosing the gauge in Eq. ( A15 1 is that it reduces to the flat gauge in the isotropic limit where -0 = -0. 



Appendix B: Integrating out non-Dynamical Fields 



In this Appendix we present the detail analysis of integrating out the non-dynamical fields SAq, /3, A and B in terms 
of the dynamical fields 6p,'^,6Ai and M. The second order action is given in Eq. (30 1. Correspondingly, the second 
order action for the scalar perturbations in Fourier space is 
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drjcTk 



bb'kl{A*(3 + An + y + ^-^)kl{A*B + AB*) + "^klkl{j*A + jA*) - a'b'V{p,)\A\' 



~b'plAl\Af - "^klkliP*B + PB*) + "^klkli^*^ + 7/3*) + ^fc^fc^(7*/3' + iP'l + ^b^plAlPM' 



"-klklW - \klkl[B*^' + By*) + \{\- "^)klkl{j*B + 75*) + jklkl\B\' + jklklljf 



~b^ AlplklklH^ 
62 



2l2 



ab I 



4 ' b 
62 h" 



C-T - -)fc^fc?l7p + ^l-Jpf - ^-^p'oiA*Sp' + ASp'*) 



4 ' 6 



p'okii(^*Sp + f3Sp*) - -p'oeJB*Sp + Bdp*) - -ki\6p\ 



^-kl\5p? 



-pI\5A,\ 



+ik,^plA,{P*5A^ - /3MS) - ^pI\5A^\'' - ^Al\6p\' - e262p„yl,(VMi + 5pSAl) 



+ik^kl—plA^{-lM* - 7*Af) 



e a 



e262 



6^62 



plklW]-" - —pIAM*5A, + AMD - ^P^Al{A*5p + A<5p*) 



2 "^'w 

62/2 



2 j-2 



{6A'*5Aq - 6A[6Al) 



-ik. 



6V: 

' 2a? 



■A;(A*Mo - A5AI) - ik^kf^^{-fM'* - -f*M') + ik^kf-^A'^i^SA* - 7*Mo) 



2a2 

fc2 

^ 2a 



2a 



2 6f2 

Aifc2(BMt + B*SAi) + ik^^A'kl(B*M - BM* 
" 2a 



f2 f2 



kl\6Ao\ 



P 



kl{M'*5Ao + M'SA*) - ^kl\SAi\^ - l^klk^M]'' + ik,'—k'y{5AlM - 5A^M*) 



2 



^ A'{5A'*5p + 5A\5p*) + ik 



P A' / X A* X „ X A X„*\ b f f ,p .,2/ A* 



^^II^A:,{SA*Sp - SAoSp*) - ^!/|i^A^^(A*<5p + A6p*) 



2a2 



P /l'2|r |2 I b?ff,pp 4/2|c |2 



a'b' ,2 a2 



2a2 

'-Vp{5pA* + Sp*A) 



(Bl) 



We have to integrate out the non-dynamical variables {6Aq,I3,A,B} from the action Eq. ( |B1[ ). The analysis are 
simple but tedious. To outline the analysis, here we demonstrate how to integrate out /3. The action expanded in 
powers of /3 is 



C = ci/3/?* + C2/3* + c;/3 + 



(B2) 



in which the dots indicates the rest of the action containing the dynamical fields {5p, 5Ai, M, 7} and {(5^0, A, B} and 
C1.2 are functions which can be read off from the action Eq. (Bl) 



ci = "^klkl + ^b'p'klAl 



(B3) 



and 



C2 = 6 6'fc2A - ^fc2^2 B + "^klkl 7 - klkl 



ab\ ' 



62 e262 

—p'klSp + ikj.—^p'^AxSAo . 



4 -x-y- . 2 

Varying the action with respect to /3* yields /3 = —C2/C1. Plugging this into the action yields 

|2 



c = - 



\C2\ 
Cl 



(B4) 



(B5) 



Following the same steps to integrate out SAq^A and B we can write the dynamical action as 



L 



(2) 



MM 



LaiAi + Lp^ + LpM + LpAi + L^M + LjAi + LmAi 



(B6) 



in which the dots indicate the rest of the action coming from the dynamical fields {5p^5Ai,A'I,^}. Here we have 
defined 
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Lpp - 



1^ 
As 4 



Po 



AiaAi 4 



62 



2 7,2 



Ai 4 ^0 A2 



lAsI^ lAgp |Ai4p 



As 



Al3 
,2j>2 a2 



(5p 



Agfe^ . _ A6Ai4 

Ag 2 AgAis 2 



IA--I2 



4''^''^ 16 A2 ^ 16 AsAi ^ ^ Ai3 



Po 

9 



{5p6p*) 



(B7) 



r, ^^^xPO^^ky + f^k^ky 



b^b_ 

4 ^ 6 
A11A12 

Al3 



« ^'.2,^2 ^ 1 -.4 4'2,2,4 _ (A4) Ay 



A« A 
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(77*) 



9 



2,2 a^-^7A3 2>2 
4 A8A2 ^ ^ 



(B8) 



.hi-L^kt- 



f |A3|^ 



j^4 _ |Al7| 



M 



2aAi3 



4 AgAi y Ai3 
fA^k^kU {i\nM*M' - i\\^M'*M^ 



2a2- 



+ 



2As 



4a4A2 
'-B^p'A, 



(Aio)2 (Ai6)2 



^13 



Ml 



^^Aio 2 2 ' AiflA 



15^16 
Al3 



I I *^ 2 2 7,2 J 1,2 1,2 



4a2Ai3 



„2 2 ; 2 

-e p 



4 4^2 (Ais)^ 

2 4A8 ^ ^ Ai3 



M 
(B9) 
2 

SAi 
(BIO) 



l2a3Ai A2 Ag Ai3 ^ V ^ V 4 A2 ^ 



('2; 2; 2 A4A5 AyAg A12A 



ab Ar 



ab A3A9 2 2 
4 A2Ag^^^^ 



AiiA 



14 



A13 

6' AgAn 
2 AgAis 



62A7 '_6^A6Ai2 

2 Ag 2 A8A13 



'a63 A3 



p)(<5pV+c.c.) + (^^pfc^fc^ 



p ) (Spj* + c.c.) 



(Bll) 



C.C.) - (^) (a..,*m' + C.C.) + C^^j^/A^pk^kl) (^SpM* + c.c.) + {^jXp) 
(A3VM'* + C.C.) - ) (a W*M' + C.C.) 



(B12) 



L,^, = (-eVpA. + ^eV^x-^)(WI+c.c.) + (^A;//, 



A14A1 



A9A10 A14A16 



As 

'^^Aio 
2 Ag 



A 



13 



Ai3 

Sp*6A[ + c.c. ) + 



2o2 A2 



4Ag 



kxf 



c.c. 



2 AgAi3 



a6A3 



- {7;^^fklkt](XslM*+c.c. 



•■8X2 A2Ag 



(B13) 



(B14) 
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^''^'"^(j'SA^+cc) (B15) 



A 



and 



= ( - Ifkkl + ^fj^XH) {ham- + c.c.) - ( j^/'fc*:?) (iM'.-M' + c.c.) 

(A3M1M'* + C.C.) - ^ (AJyMiM'* + C.C.) + ^^J^^xkxkl (i5A-,M* + c.c.) (B16) 
The parameters Aj and Aj which are introduced after integrating out the non-dynamical fields are defined via 

Ai =^fcV^ + y&V (B17) 

A2 =\klkl^'-l^fp'^Al (B18) 

A3 =hh'kl-^^fp-A^A'^ (B19) 

X. =^(^^>^.?.^/V...; (B20) 

A5 =-|-p'fc^ + |0//,pP^A.A; (B21) 

Ai =^klkl + ^klp'Al (B22) 



A2 =^^^ + ^A^^^ (B23) 



A3 =-lg^/2yl^ + i^fc3^2^^ ^B24) 

A4 =-^^^ + ^(-+3^j^^fcyA. (B25) 

As ='4kx.ff,,K-'-^klp'p'A, (B26) 

a6,9,2 abH' f a' fo'\,4,2 A5A4 



As 


= -aH'^V - 


626'2 IA3P e262 
Ai '"^ A2 2 ^ 


A9 






Aio 


~ 2a? ^ ^ 


*^'^3 .2, 
2a2A2^ " 


All 


= ~^klky - 


a'&',4,4 «'>A3A5 2 2 
16X2 ^^ '' 4A2AS ^ ^ 
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52 /a' 


L2.2 1 «&A4 2,2 A7AS 

6 J "^''y + 4A2 ^ As 



„2i,2 



(B29) 

pAl (B30) 

(B31) 
(B32) 
(B33) 
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7 9 T 9 1 \ 1 9 

"4''-''^ I6A2 As 


(B34) 


Ai4 


_ , «^A5 2,2 AgA^ 
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(B35) 
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J 9 1 9 \ 


(B36) 


Al6 


Ag lDaAiA2 


(B37) 


Ai7 


- 2A2As^ ' 16AiA2^ ^-^^-^ 


(B38) 



As mentioned in the main text, amongst A,; and Ai, the parameter Ai plays the key roles which determines the form 
of other A^ and A^. One can check that Ai originates from integrating out 5Aq. 



Appendix C: Second order Slow roll action 

In this Appendix we calculate the second order action for the canonical fields in the slow roll approximation. 



Following the discussions in Section III we divide the dynamical action into two different regions depending on 
whether the charge e is important or not. In the first phase, 77 < rjc, the charge effect is sub-dominant while in the 
second phase, rj > 77^, its effect is dominant and the longitudinal mode, as we will introduce it in Eq. (C12), has the 



same contribution as the transverse mode. In the following, first we write the action in the first phase and then we 
go to the second phase. 



a. First Phase 



Our goal is to write down the action in terms of the free Lagrangians plus the interaction terms. Using the slow-roll 
approximation given in Eq. ( 35 1 we have 



Where 
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(CIO) 



(Cll) 



Now looking at Eq. (Cll I we can easily see that this term is not as small as the other interaction terms. Actually this 



term seems to be of the same order as our free field action. This means that M and SAi are not the physical fields. 
One should consider a rotation in {M,6Ai} space such that all of the interaction terms become small compared to 
the free field action. One can easily check that the following two new fields Di and D2 work for us in the sense that 
they do not mix with each other and all of the interaction terms would be small: 



Di = 5Ai - ik cos 9M 

D2 = cos 96 Ai + ik sin^ 9M . 



(C12) 
(C13) 



One can check that Di represents the transverse polarization while D2 is for the longitudinal polarization of the gauge 
field perturbations SA^. Now the different parts of the action can be rewritten in terms of these two new fields as 



,2 f2 



LsAiSAi + LmSAi 

LspM + LgpSAi 
L^M + L^SAi 
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(C16) 



where in Eq. (C15) and Eq. (C16) only the leading terms have been written. It can easily be seen that since the 



charge effect is not important in this phase, the longitudinal mode is sub-leading. Now we can define the canonical 
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variable as, 



5p — bSp 



7 



-/sin 9Di 



2V2Xk y l~I 



(C17) 
(C18) 

(C19) 

(C20) 



One can write the action in terms of the canonical variables. Due to technical reasons we only express the free 
Lagrangians in terms of the canonical variables while the interaction terms may be expressed in terms of the old 
variables 



L. 

L 

LD2D2 
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DiDi 



1, 
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(1 - 2 sin^ 61) 



= M 
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1 - / 
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1 
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\D2 



(C21) 

(C22) 

(C23) 

(C24) 
(C25) 



These are the final results for the action in the first phase used in Eqs. (37 1- (40) 



b. Second Phase 



In this part we write the leading order action in the second phase fro m which we can read off the canonical fields. 
Considering the dominant effects of e in Ai as mentioned in Section III yields 



Lpp 

LdiDi 

LD2D2 

LpDi 

LpD2 



&2 

2 

2a2 



1\ 



IeH\\ 



Af4 ) \prf 
16A2/2 



,f8inH]\Di\ 



16A2/2 



1\ 



6Ism^ef[Sp*D^ +C.C 
QI cone f {5 p*D2 + c.c 



JlJ ^ 

Now from the above equations, we can find the canonical variables in the second phase as, 

5pk = bSpk 

Dik = - sinOfDik 
a 

Ihk = -fD2k 
a 



(C26) 
(C27) 
(C28) 
(C29) 
(C30) 

(C31) 
(C32) 

(C33) 



as used in Eq. (59). 
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